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ABSTRACT 

A tomographic method is described to quantify the three-dimensional power-spectrum of the ionospheric 
electron-density fluctuations based on radio-interferometric observations by a two-dimensional planar array. 
The method is valid to first-order Born approximation and might be applicable to correct observed visibilities 
for phase variations due to the imprint of the full three-dimensional ionosphere. It is shown that not the iono- 
spheric electron density distribution is the primary structure to model in interferometry, but its autocorrelation 
function or equivalent its power-spectrum. An exact mathematical expression is derived that provides the three 
dimensional power-spectrum of the ionospheric electron-density fluctuations directly from a rescaled scattered 
intensity field and an incident intensity field convolved with a complex unit phasor that depends on the w-term 
and is defined on the full sky pupil plane. In the limit of a small field of view, the method reduces to the single 
phase screen approximation. Tomographic self-calibration can become important in high-dynamic range ob- 
servations at low radio frequencies with wide-field antenna interferometers, because a three-dimensional iono- 
sphere causes a spatially varying convolution of the sky, whereas a single phase screen results in a spatially 
invariant convolution. A thick ionosphere can therefore not be approximated by a single phase screen without 
introducing errors in the calibration process. By applying a Radon projection and the Fourier projection-slice 
theorem, it is shown that the phase-screen approach in three dimensions is identical to the tomographic method. 
Finally we suggest that residual speckle can cause a diffuse intensity halo around sources, due to uncorrectable 
ionospheric phase fluctuations in the short integrations, which could pose a fundamental limit on the dynamic 
range in long-integration images. 
Subject headings: instrumentation: interferometers — techniques: interferometric — methods: analytical 



1. INTRODUCTION 

Over the last years, considerable effort has gone into the 
construction of large wide-field low-frequency radio interfer- 
ometers. One of many challenges that these new instruments 
must overcome is that of phase and amplitude fluctuations 
due to the ionosphere (Hewish 1951, 1952). This is partic- 
ularly important at meter- wave legths whe re these instruments 
operate (see e.g. lCohen & Rottgering l20"09) and the refractive 
index of the ionized medium can be particularly larg43- More- 
over, low-frequency radio telescopes often have wide fields of 
view (i.e. up to tens of degrees), making them not only sus- 
ceptible to the combined electric field from large parts of the 
sky, but also to its modifications by the full three-dimensional 
electron-density structure of the ionosphere. 

For high-frequency interferometers - having dish-antennae 
that probe a small area of the sky (a few degrees or less) - 
the three-dimensional ionosphere can be collapsed along the 
line-of-sight and be well approximated by a two-dimensional 
thin phase-scre en, without loosing substantial information 
(ISalpeterlll967h . This approximation, however, does not hold 
if the field of view becomes large and widely separated an- 
tennae observe different sources, under large angles, through 
similar parts of the ionosphere (e.g. Lonsdale 2005). In 
the case of wide-field low-frequency interferometers, a full 
three dimensional model of the ionosphere is required for 
self-consistent ionospheric calibration. Whereas considerable 
work has recently gone in to developing (mult i-layer) two- 
dimensional wide-field ionosph eric models (e.g. llntemaetal.1 
120091: iMateiek & Moralesl l2009l) . no physically intuitive and 
self-consistent description seems readily available for full- 

' The ratio of the plasma over the observing frequency at 150 MHz is 
~0.03; hence intensity scattering is at the level of ~0.001 and scales with A~. 



scale three-dimensional wide-field modelling of the iono- 
sphere. 

In this paper a step towards such a complete physical 
description is presented based on the tomographic theory 
of electric fi eld scattering by weakly inhomogeneous media 
(IWolflll969l) . as applied to the ionosphere. The main result 
is that the ionosphere, over a wide field of view, acts as a 
scatterer with a spatially varying point spread function. The 
instantaneous two-dimensional point spread function at a po- 
sition s (a directional unit vector) around a point source at 
So is identical to the instantaneous three-dimensional power- 
spectrurrQ of the ionospheric electron-density sampled from 
points (s - So) on th e surface of an "Ewald sphere of re- 
flection" (E waldlfT969 ). i.e. a spherical surface with unit ra- 
dius centered on the direction of the source. Having many 
point or very compact sources over a larger field of view al- 
lows one to sample the three-dimensional Fourier structure 
of the ionosphere at many distinct points. These points can 
be used to constrain a three dimensional model of the power- 
spectrum of the ionosphere. This model can subsequently be 
used to "deconvolve" the image, also at points where there 
are no strong sources, and remove the dominant phase errors 
due the three dimensional ionosphere. We term this "tomo- 
graphic self-calibration", because it involves a three dimen- 
sional ionosphere and not a single phase screen. This cal- 
ibration can not be done in th e classical way through self- 
calibration (see e.g. Pearson & Readhead 1984) because the 

^ Throughout the paper, if not mentioned otherwise, the power- spectrum is 
assumed to be the exact squared amplitude of the ionospheric electron density 
waves as function of their three dimensional wave vectors at a given moment 
in time, and not the expectation value of this squared amplitude over an en- 
semble average IBourgoiJ1981D . In general the former is a Gaussian random 
realization of the latter. 
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scattering point spread function is not spatially invariant as 
is implicit in that method. It either re quires the solution of 
a matrix (m easurement) equation (e.g. ' Hamaker et all 1 19961: 
iMatejek &" Morales 2009; Liu et al. 20 IC^, which is computa- 
tionally more demanding than traditional self-calibration or, 
as we will show, it can be written as a specific three dimen- 
sional Fourier transform. 

The main goal of this paper, however, is not to introduce a 
specific solution scheme or algorithm to these equations, but 
to provide more useful guidance to future methods of self- 
consistent physics-motivated three dimensional modeling of 
the ionosphere and its use in calibration. The second goal is 
to describe several effects of a three dimensional ionosphere 
on interferometric measurements at low frequencies that go 
beyond thin phase-screen models, which have thus far been 
very successful, but are demonstratively incorrect for an ex- 
tended three dimensional ionosphere. 

The outline of the paper is as follows: In Section 2, the the- 
ory of scattering of a plane-wave electric field is restated in 
terms more familiar to radio interferometry and extended to 
the case of multiple point sources. In Section 3, the cross- 
correlation of the scattered electric field is determined (i.e. 
the visibility function) and several effects of the ionosphere 
on imaging are described (e.g. speckle). In Section 4, the first 
order effect of the thickness of the ionosphere is analyzed. In 
Section 5, the description is extended from multiple points 
sources to a continuous intensity field and a mathematical ex- 
pression is derived that allows one to build a three dimen- 
sional model of the power spectrum of the ionosphere from 
information obtained only in the interferometer plane (i.e. a 
holographic principle). In Section 6, it is shown how the 
tomographic method, connects to an extension of the phase- 
screen approach to three dimensions using a Radon transfor- 
mation and applying the Fourier projection-slice theorem. In 
Section 7, we summarize our results and give conclusions. 

2. WEAK SCATTERING BY THE IONOSPHERE 

In this section, the basic theory of electric-field tomogra- 
phy of a weakly scattering media is restated as first intro- 
duced by (Wolf 1969). The notation is adapted to that typ- 
ica Uy used in radio a stronomy and radio interferometry (see 
e.g. [Thompson et £1112001) . Under the assumption that the re- 
fractive index of the ionosphere varies slowly over a single 
photon wavelength, and the ionosphere is "frozen" over the 
time scale the radiation passes through it, one can write the 
electric field equation in a dielectric medium as 



V^E„(r) + [^X(r)]Ev(r) = 0, 



(1) 



where ny(r) is the refractive index of the medium at a posi- 
tion r and at frequency v, and k - 2n/A is the usual wave 
number for a wavelength A - c/v (iBorn & WoH)ll999h . In 
a plasma with electron density rteir) the refractive index is 
given by n^ = 1 - [nee^/iv^nieeo)] = 1 - (vp/v)^, where Vp is 
the plasma frequency, being typically ~5 MHz for the iono- 
sphere. We drop the explicit frequency dependence of the 
electric field, but it should implicitly be assumed. 

Because the three components of the electric field are inde- 
pendent, each component satisfies the same solution and we 
can simply use the scalar E to describe the electric field. If 
the refractive index of the ionosphere is near unity, which is 
often the case in radio astronomy (i.e. v » Vp), the equation 
for the electric field can be conveniently rewritten as 



where <l)(r) = k^[n^ - 1]/4jt is called the scattering potential. 
In case of n »; 1, the scattering potential strength is close to 
zero and scattering is wealfl. When <l'(r) - everywhere, the 
equation reduces to the Helmholtz equation for a plane wave. 
The s olution of this equation can be obtained dBorn & WoJB 
119991) through Green's functions and leads to an implicit inte- 
gral equation 



E(r) = E'^ir) + f <D(r')£(r')- 

Jv |r 



7(1) ( 



ik\r-r'\ 



r'l 



-dh'. 



(3) 



where ^(r) = £''^''(r) + E^'^\r). The first term is the incident 
(plane) wave and the second term the scattered wave. The lat- 
ter is equal to the integral term above and carried out over the 
entire volume V in which n i^ I. The last factor in the inte- 
grand indicates that £'*''(r) is a spherically outgoing wave, as- 
suming the extent of the scattering potential is small compared 
to the distance between observer and scattering medium. 

2.1. Weak Scattering of Multiple Point Sources 

We now extend the single plane-wave description of lWoHl 
(119691) to an incident electric field that results from the sum 
of N point sources that satisfy the solution of the Helmholtz 
equation in free space. Later in the paper we further extend 
this to a continuous intensity field. We also express all phys- 
ical distances in units of A, i.e. u = r/A = (u, v, w) in the re- 
mainder of the paper. Its Fourier equivalent is s = (s„, s,,, s„)- 
The incident electric field in this case becomes 



£('\u) = 2 V5^^< 



2;r/So,„-u 



(4) 



(V^ + r ) E(r) = -47T 0(r) E(r), 



(2) 



where S„ is the flux-density of point source n - I . . .N and 
So,« are unit vectors that point in the directions of the point 
sources. These point sources are (for now) assumed to dom- 
inate the electric field and can be compared to the phase cal- 
ibrators in radio interferometry. In the weak scattering ap- 
proximation, the scattered wave has a relatively low amplitude 
compared to the incident wave. We can then replace ^(u) with 
^^''(u) to first-order (Born) approximation, finding 

Zf „2m|u-u'| 

V^; <i.(u')e2-^«»-'^ — -^'u', (5) 
Jv lu-u'l 

where the subscript indicates that the scattered field is that 
to first order Born approximation and <t>(u) = [n^ - 1]. We 
note here that the second order can be obtained by substituting 
^''^(u) + iij (u) in to the original integral equation to obtain 
a solution to second order, etc. This iterative scheme only 
works for weak scattering, where the values of 'I>(u) do not 
exceed unity. 

If we us e the fact th at the spherical outgoing wave can be 
written as (IWevll919h 

|u - u'l J_oo J_„ s„ 

(6) 

^ We note here that even though scattering can be "weak" in terms of a 
small deviation of the refractive index from unity everywhere, the integral 
over the line of sight can still lead to substantial phase changes in the electric 
field, leading to both dift'ractive and refractive effects and strong scintillation 
(e.g. Naravan 1992). In the current discussion "weak scattering" means both 
a refractive index very close to unity and phase fluctuations much less than 
unity. The latter is the nominal mode of the ionosphere and data obtained 
with interferometers during strong scintillation (e.g. during solar bursts or 
sun rise/set) occur only rarely and are often discarded in analyzes. 
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with si. - I - si + si for si + si < I or si - -(si + si - I) 
otherwise. Complex values of s„, lead to exponentially de- 
caying (evanescent) electric fields that are typically not mea- 
surable far from the scatterer The other (homogeneous) 
waves are those measured by a distant observer If one fur- 
ther uses the Fourier transform of the scattering potential, 
3)(s) = JjJ <I)(u)e"^'^'* "c/^u, one can write the scattered elec- 
tric field as 

Efiu) ^iy^/s~„ff — 0(s - So,„) e^^'^-^ds^ds,, (7) 
JJ s,, 

where we assume a geometry where w = is the ground- 
plane below the ionosphere where n = 1 , and that w > is in 
the direction of the zenith or the phase reference center (see 
below). The interferometer is placed in a plane defined at a 
constant Want - Zant/^- Typically one can assume Want - 0. 

Thence, one finds a relation between the Fourier transform 
of the observed electric field in the plane of the interferometer 
at Want and the Fourier transform of the scattering potential 



E^'Hsu, .„) - -e+2-....n„„ -^ ^3,(, _ ,^ .)^ 



(8) 



with E^'\su,s,) = ]J"£f(M,v,Want)e+^'''<'""+''"^rfMt/v. This 
can be regarded as the Fourier transform of a two-dimensional 
slice through a three-dimensional scattered electric field. In 
this paper we do not treat the case of an interferometer with 
varying Want- A planar array is an reasonable assumption for 
relatively compact (i.e. km-scale) interferometers, but breaks 
down on large scales where the curvature of th e Earth can not 
be neglected (see e.g. ICarozzi & WoanI [20091) . For a planar 
array, however, the w-term due to the array can be neglected 
for small integration times (i.e. instantaneous sampling of the 
electric field in a plane), in contrast to visibilities from very 
different time frames where the array has rotated over a sub- 
stantial angle compared to the phase center (only a linear east- 
west array does not suffer from the w-term). 

The physical interpretation of Eqn.® is the following: Ev- 
ery point of the two-dimensional Fourier transform of the 
scattered electric field in the plane of an interferometer probes 
a single three-dimensional mode of the scattering potential 
(i.e. the scattering medium) for a single point source. In the 
presence of A^ point sources, all in different directions, every 
point of the two-diminsional Fourier transform of the scat- 
tered electric field in the plane of an interferometer probes the 
sum of A^ independent three-dimensional modes of the scat- 
tering potential. In Section 5 we show how to unravel this 
information. 

3. THE EFFECTS OF WEAK IONOSPHERIC SCATTERING ON 
INTERFEROMETRIC IMAGES 

In radio interferometry one does not analyze the electric 
field itself. In that case, Eqn.® would directly yield the three- 
dimensional structure of the ionosphere (per integration time) 
because the phase information of the Fourier transform of the 
electron density of the ionosphere is fully retained in the phase 
information of the scattered electric field. In reality, only the 
cross-correlations of the electric field, measured at different 
antennae pairs, are stored (i.e. the complex visibilities) and 
the phase information of the ionospheric density fluctuations 
is lost. In the following, we assume that the total electric field 
from the entire sky (i.e. the antenna sensitivity is direction- 
ally independent) is measured over the infinite interferometer 
plane with w - Want- 



Visibilities are sampled from the cross-correlation of the 
electric field ^(u) = ^'''(u) + E^^\u) with its complex conju- 
gate, i.e. V(b) = {E{u)E*(u + b))t with b being the baseline 
between two points (antennae) in plane of the interferome- 
ter. The averaging is assumed to be over time. The Fourier 
transform of the visibilities forms the incident intensity from 
the sky, as follows from the van Cittert-Zernike theorem (e.g. 
ICarozzi & Woanl2009) . The same intensity is also the product 
of the Fourier transform of the electric field with its complex 
conjugate. A bit of algebra shows that the cross-correlation 
between the incident and scattered fields depends on the imag- 
inary part of the zero-mode, 3)(0), of the ionosphere, and con- 
sequently is equal to zero. The multiplication of the Fourier 
transform of the scattered electric field with its complex con- 
jugate therefore provides the complete scattered intensity 



6I^%„ s,) = <£*-'>*(i,„ s,)E^-%,, s,.))u 



(9) 



where the dependence on Want disappears. Using Eqn.®, we 
find the following result 



<i(^«,^v.) 



— J] J] ^JS nS ,n^*{S-S0,n)^iS-S0,,„). (10) 



This equation is exact for phase-coherent point sources to first 
order Born approximation. Howev er, the sky is an incoherent 
emitter (see lMandel & WoI3ll965L for an exposure on the co- 
herence properties of electric fields). Hence, the cross-terms 
with n i^ m depend on the electric field coming from incoher- 
ent point sources and vanish, such that we are left with 



5I^'\su, s, 



>')=iZ^"l^(^- 



so,«)r 



(11) 



where we dropped the subscript. This equation forms the ba- 
sis for further discussions in the paper. The above equation is 
only correct for an interferometer and an electric field mea- 
sured in a plane. In three dimensions, one would no longer be 
able to use simple Fourier transforms (see below), because s„ 
depends explicitly on Su and s^. 

To understand the physical interpretation of the above equa- 
tion, one might suppose a point source in the zenith (or equiv- 
alently in the phase center) emitting a plane wave in the ab- 
sence of the ionosphere. Because the phase of the electric 
field is the same at each antenna (by construction), its Fourier 
transform yields a complex delta function in the zenith with 
a time-varying phase. Multiplied with its complex conjugate, 
this recovers the point source intensity. If a two-dimensional 
thin phase-screen is placed in between the source and the ar- 
ray, exhibiting a single wave-mode in electron density perpen- 
dicular to the zenith or phase reference center direction, then 
part of the electric field amplitude will be modulated such that 
its phases show to first order the imprint of this ionospheric 
wave-mode (see RatcUffe.1956 . for a wonderful description). 
The modulated phase (i.e. a single wave over the array) can 
be interpreted as being identical in the weak scattering limit 
to the modulated phase of a point source offset from the zenith 
in the direction of the ionospheric wave-vector by a distance 
set by the phase-frequency over the array. 

Hence, a single ionospheric wave-mode scatters a fraction 
of a point source intensity into an additional point source (a 
speckle), offset by the projected phase frequency onto the ar- 
ray and with an intensity proportional to the amplitude of the 
wave-mode squared. The sum of all speckles create a halo 
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of scattered emission around the point source, when not cor- 
rected for through phase calibration. 

3.1. Speckle and Speckle Noise 

Because the phase fluctuations of the electric field are mea- 
sured in the plane of the array and are typically drawn from 
a Gaussian random field realization from some (ensemble av- 
erage) power spectrum, taking the Fourier tranform of this 
field and multiplying it with its complex conjugate yields that 
each point source exhibit s a diffuse "h alo" of scattered inten- 
sity (e.g. S alpeter 1967: Cronynlll972h . Its Fourier transform 
is related to the usual phase structure function (D^; see be- 
low). This pattern is referred to as "speckle" in optical (laser) 
interferometry. We define the instantaneous ionospheric scat- 
tering PSF (ISP hereafter) as the sum of a delta function plus 
its scattered speckle pattern, re-normalized to a flux-density 
of unity such that flux is conserved in a convolution process. 
The ISP i s related to the Fourier transform of optical transfer 
function JGoodmanl flQSS). determined by phase fluctuations 
induced by the ionosphere. Eqn.dTTTl shows that the point- 
source intensity is convolved with an ISP that reflects a curved 
surface in the Fourier transform of the instantaneous iono- 
spheric electron-density fluctuations. Disregarding a geomet- 
rically determined distortion that depends only on the angle of 
the observed point on the sky away from the zenith, a three- 
dimensional ionosphere causes a spatially varying convolu- 
tion. In contrast, a thin two-dimensional ionospheric screen 
causes a spatially invariant convolution. This sets our analysis 
apart from most studies up to the pre sent that have fo cussed on 
scattering by ionized media (see e.g. lNaravanI 19921 for an ex- 
cellent review) where wide-field effects and extended (thick) 
screens can be neglected in nearly all circumstances. 

For ionospheric electron-density fluctuations set by Kol- 
mogorov turbulence, one expects that |€)(As)p oc lAsr'^ with 
11/3 over a scal e of meters to tens of kilometers (e.g. 
[Thompson et al.ll2001h . This implies that in the weak scatter- 
ing regime, point sources on average exhibit a speckle pattern 
which rapidly decrease s in intensity with distance from the 
source (Goodman 1985). The speckle pattern, at a given mo- 
ment in time, is a Gaussian random realization from a power- 
spectrum with expectation value |0(As)p. Scattered flux is 
spread over scales from arcseconds to degrees, correspond- 
ing to tens of km to tens of meter-scale ionospheric electron 
density fluctuations, respectively. It is interesting to note that 
ionospheric modes of a gives physical scale give rise to image 
distortions that are visible only on baselines equal or larger to 
that same physical scale. The properties of large isoplanatic 
patches over a given scales are therefore easier to determine 
if one can observe the sky with an array that exceeds that 
scale (see also e.g. Jacobson & Erickson 1992). Long base- 
lines thus significantly help in calibrating shorter baselines, 
even if the science of interest is obtained on the short base- 
Hnefl 

The presence of a speckle pattern can have consequences. 
Once the ISP reaches an intensity level comparable to the 
average intensity of the surrounding sources, the confusion 
noise, or the noise itself, one can no longer distinguish it from 
these contaminants and calibration is typically limited to a 
maximum ^-mode. This could place a limit on the calibrat- 
ibility of the instrument, especially at very low frequencies 

^ To correct for all modes in the field of view (6fov). the longest baseline 
should be fc > 25 km X [//ion/(300 km)] ■ [eF0v/(5 deg)], where //!„„ is the 
typical height of the ionoshere. 
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near the ionospheric cutoff", because the signal to noise ratio 
of an image snapshot, taken over the maximum time-scale of 
substantial ionospheric changes, might not be large enough to 
correct for low signal to noise speckles further away from the 
sources. This leaves residual speckles which are a function of 
wavelength (scaling with A'- in strength and with A in scale). 

A detailed calculation of the residual speckle intensity and 
its associated noise goes goes beyond the scope of the cur- 
rent paper. We note that these speckles do not average away 
over time but ultimately form a halo (i.e. "seeing") around 
each source. Whereas a bright source might be very useful 
in determining the ionospheric scattering, it might also leave 
a residual speckle pattern in the image due to uncalibrated 
phase ffuctuations from ionospheric scale below the smallest 
baseline. It could therefore be advantageous to place bright 
sources near the half -power of the field of view such that they 
still probe the same ionosphere to first order, but that their 
residual speckle patterns do not contaminate the field of inter- 
est. 

Speckle noise (i.e. the expected standard deviation from the 
expectation value of the speckle intensity at a given point), 
however, does average away. To see this, we note that a con- 
volution of the intensity implies a multiplication of the visibil- 
ities with the Fourier transform of the ISP. In the case of weak 
scattering and assuming a thin phase screen and small field 
of view, this multiplication function to first order is I - Dd,{0) 
where the latter term is the phase structure function. ' Cronvnl 
(1972) showed that the variance in the scattered part of the 
intensity is set by the amplitude of D^{6) is to first order equal 
to the its expectation value. Hence speckle noise around point 
sources average away as one over the square root of the num- 
ber of independent realizations of the ionosphere. Obviously 
calibration can remove much of the speckle pattern, but not 
that caused by the structure below the smallest and above the 
longest baselines or below the noise level, respectively. This 
leftover halo of speckle can be compared to that in optical (af- 
ter AO correction ) in the search of fa int companions around 
nearby stars (e.g. iRacine et alj|1999l) . Speckle and speckle 
noise ultimately set the detection threshold and a similar ef- 
fect can play a role in very high dynamic range imaging at 
low radio frequencies. We finally note that a convolution im- 
plies a suppression of visibilities on long baselines. In partic- 
ular, if no calibration is applied for ionospheric effects, phases 
on long baselines beco mes less correlated and average out 
stronger jBramlevI 19541) and a "seeing" halo develops around 
point sources. 

4. VARIATION IN THE IONOSPHERIC SCATTERESIG PSF OVER THE 
FIELD OF VIEW 

In the case of a single source in the phase-center and narrow 
field of view over which 5, changes little from zero, the ISP 
around a point source (assuming a perfectly sampled uv -plane 
or electric field) represents the instantaneous power-spectrum 
of the electron density ffuctuations on a two-dimensional sur- 
face in three-dimensional Fourier space. In the case of a nar- 
row FOV and a compact speckle pattern with ,Svj, — > 0, this 
pattern is nearly equal to the instantaneous two-dimensional 
power spectrum of the medium integrated along the line-of- 
sight, which has thus far motivated the use of a phase-screen 
or thin ionospheric model (e.g. Intema et al. 2009). 

To assess how the ISP is modified when displac ed from the 
phase reference center (see lThompson et ani200lL usually the 
center of the antenna beam), we assume that the phase refer- 
ence center is in the zenith, such that the w-axis points exactly 



Ionospheric Power-Spectrum Tomography in Radio Interferometry 



upward and the Mv-plane is perpendicular to it. The planar ar- 
ray lies at constant w as assumed thus far. It can be shown, 
however, that for any planar array, the coordinate system can 
be rotated to any phase reference center By correcting the 
phases of each of the antennas such that the fringe-rate is zero 
for a chosen phase reference center, one effectively re-orients 
the planar array such that it acts as being perpendicular to 
the line connecting the array and the phase-reference center. 
Hence, the analysis performed in the coordinate system with 
the phase reference center in the zenith is in fact valid for any 
chosen phase-reference center for planar arrays by defining 
the w-axis in that direction. 

We now assume two point sources, one in the phase refer- 
ence center at Szenith - (0, 0, 1) and one slightly offset from the 
zenith at a unit vector Ssi-c, such that |Ssi-c - SzenithI **: 1- One 
can then Taylor expand the ISP to first order To do this prop- 
erly, however, one needs to compare two points that are offset 
from the two sources (on in the zenith and one offset from it) 
by an identical two-dimensional vector 5S2D = {5su,5s^). In 
that case, one compares the ratio of their aligned speckle pat- 
terns in the (i„, i,.) plane. Furthermore, because the first two 
coordinates of s - Ssrc = {5%2-d,5sk - Ss'^'^ = (1 - |5s2dP)'^^ - 
(1 - I^SjQp)'''^) are the same by construction for both sources, 
any difference occurs because of the dependence of the ISP on 
6s^^!^ and hence one only needs to Taylor expand Eqn.lfTTTi for 
a single point source with respect to changes in Ss^^, which is 
the offset of the source from the phase-center 

After a bit of algebra one finds that the fractional intensity 
difference between the offset ISP divided by the one in the 
phase center is 



f{s) 



/;^'(5s2d)^ 



2 8<^^'^ ■ 



6S2D 



1 - |5S2dP 



1 



Vi - |5s2dP dm)? 



2|(D(t)p 



du 



(12) 



with the derivative of the power-spectrum being evaluated at 

the point t - (8^20, V1-|^S2dP - 1). The length of 5s2d does 
not necessarily have to be much smaller than unity because so 
far we only expanded in Ss^^. If we further assume that also 
I^S2dI ^ 1 then the equation further simplifies to 



fU) 



/;'(<JS2d)^2 5s^55-5S2D 



1 



51n mi)? 



2du 



(13) 



The ISP is therefore no longer spatially invariant but starts 
to depend on the power-spectrum of the ionosphere in the w- 
direction and on the offset of the source from the phase refer- 
ence center. Whereas to first order the ISP remains identical to 
that of a thin phase-screen, the large-scale structure (i.e. with 
t„ a: 0) of the ionosphere in the w-direction becomes impor- 
tant because d\^{i)? I dt, + constant for a thick ionosphere. 
Because this additional term depends strongest on the large 
scale w-modes, we expect this term to be rather insensitive 
to relative small changes in Sjp and show slower diurnal and 
seasonal variations. 

If, instead of working with the image 5/*'*(.s'„, .?,,), we work 
with SJ'^^Hs,,, i,,) = si, dl'-^^Su, Sy) (see also Section 5.2), the 
fractional difference is 



As) 



fy'(SS2D)^-(Ss'^^-6S2D) 



(6s^^ ■ 6s2d) 



Vi - |5s2dP dm)? 



mt)? 

51n|0(t)|2 



dt. 



du 



(14) 



A physical interpretation of this is the following; when ob- 



serving away from the phase center, one becomes sensitive to 
large scale modes in the w-direction. This causes a gradient in 
the electron density over a scale that corresponds to depth of 
the ionosphere along the line of sight to the source projected 
onto the array. For small angles away from the phase center 
this projected baseline is small and interference between the 
opposite parts of the ionosphere is not symmetric. This results 
in large scale interference patterns on the plane of the array. 
These patterns correlated strongest on long baselines and thus 
show up as small scale structure in the ISP, i.e. near the point 
sources. We note that in the case of a thin ionosphere, which 
has 51n|3)(t)p/i9f,v = 0, the ISP is spatially invariant to all 
orders. 

The interpretation of Eqn. (fl4l l is that offsets, S2d, perpen- 
dicular to Sjp lead to a zero fractional difference, whereas 
parallel offsets lead to maximal differences (either positive of 
negative). This implies that to lowest order the ISP is invari- 
ant in the tangential direction when moved from the phase 
center on radial spokes and is stretched or squeezed in the 
radial direction further away from the phase-center, depend- 
ing on the gradient of the ionospheric power-spectrum in the 
vv-direction. Eqn. (fT4l) can thus be used to determine the ISP 
over a wider field of view, taking the first-order effect of the 
thickness of the ionosphere into account. By having several 
calibrators spread over the field of view, one can determine 
the nearly constant value of d\n\^{i)? j dt„. Once known, the 
ISP can be determined to lowest order for the rest of the field 
of view by simple rescaling. 

To estimate the effect of the thickness of the ionosphere let 
us assume it has a thickness d,\ in units of wavelength and 
a uniform electron density in the w-direction. The normal- 
ized power-spectrum in s„ direction has the functional form 
|0(f,i.)p - (dxln) X sinc^(ii^fH.)- Inserting this in Eqn.(fT4b and 
Taylor expanding to lowest order yields, along radial spokes, 
ff X (7t/3) dx \5&2D? \5s^^V We note that even though the 
fractional changes in the ISP are smaller near the source^ 
the absolute changes are larger because of the steep increase 

of the power spectrum, i.e. \^{6s2d)? x ff '^ l<Js2Dr"''^ \S^2d\ 
for/? =11/3 i.e. a Kolmogorov spectrum. Note that this equa- 
tion breaks down near the outer scale (i.e. whose effect shows 
up closest to the images) of the ionosphere otherwise it would 
not be bounded. 

5. THREE DIMENSIONAL IONOSPHERIC POWER-SPECTRUM 
TOMOGRAPHY 

After having derived an equation in the previous section 
that describes the scattering of point sources due to the three- 
dimensional ionosphere, here we introduce the equivalent ex- 
pression for a continuous intensity field and then discuss the 
consequences of these results for ionospheric calibration. 

5.1. Scattering of a Continuous Intensity Field 

One can simply extend the point-source equation for scat- 
tering (Eqn. lfTTI) ) to a convolution-type operation on a con- 
tinuous intensity field by relating flux-density with intensity 

5 Numerically we find /|'' ~ lO^^ y,,/105) |5s2d/10'P |(5s™/10°|. As 
example, for an ionosphere with 200 km thickness and observed at 2 meter 
wavelength, the ISP about 10' away from sources near the edge of a 10 degree 
field of view changes by ~ 1 % of the ISP maximum. This might seem small, 
but could still be considerably larger than the noise or confusion level for 
bright sources in the field. Hence, in high dynamic range and wide field of 
view imaging experiments the third dimension of the ionosphere becomes 
very important. 
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times area. In that case one readily sees that 



df'\s,„ s,) 



U!" 



(S0,u, iO,y) l<l>(s - So)| dso,udso^y. 



(15) 

We note that this is a two-dimensional convolution with a 
three dimensional kernel. This makes it more difficult to de- 
convolve using simple two-dimensional Fourier techniques. 
For a small field of view, one can simply set s^ - so^„ - and 
perform a two dimensional convolution through fast Fourier 
transform methods. For wider fields of view this can not be 
done. However, one can rewrite the equation as a three di- 
mensional convolution as follows 



61^%,,, Sy): 



iff!"i 



(so)5k(<JiJ|0(s-So)rflfso,(16) 



where Ss^, = sq.h' ~ (1 ~ *o„ ~ ^Ov-*'^^ ^"'^ ^k is a Kronecker 
delta function. Note that in this equation So should no longer 
be treated as a unit vector but still that (1 ~ ■Sq,, ~ ■'o r-* - ^■ 
The intensity I^^(so) is a cylinder of radius unity that has the 
same value as I^'^squ, io.i') for each value of s„. Writing the 
equation like this, makes it a three dimensional convolution 
which can be performed using Fourier transform techniques, 
but at the cost of more memory and computational effort. 

5.2. Determination of the Three-Dimensional Ionospheric 
Power Spectrum 

Here we derive an expression for the full ionospheric power 
spectrum in terms of the scattered intensity field and the inci- 
dent field. To do this, first we define the rescaled version of 
the scattered intensity as 



lis) 



5r>{su, s,) = s^. ■ 6f-'>(s,„ s,). 



(")( 



(17) 



and 6P\u,v) = ff 6J'^'\s,„s,)e^^''''^'""^''''''ds,ds, as its 
Fourier transform. We note that the scatte red intensity is zero 
if si + si > 1 . If we now define (see e.g. iSramek & Schwab! 
Il989t iCornwell & PerlevI 119921 for a similar approach for 
non-coplanar array© 



4d(so) = 4',^(so) ■ 5(^0,H 



v; 



■'O.u 



So,,.), 



(18) 



we obtain 



SJ''\Su, s,.) 



Iff 



4d(so)|3)(s- 



So)rfl'so. (19) 



We note that this equation can be integrated over infinity, as 
long as the intensities are zero (as they are) when |so| > 1. Us- 
ing the relation between convolution and Fourier transforms, 
we can now write this as 



6p\u,v)^f(^(u)-r{mst)(u) 



(20) 



with f^l 



(u) = fff J%(so)e--'''^«"dsQ and :r(|0(s)|2) being the 
autocorrelation function of the ionospheric scattering func- 
tion. This remarkable equation shows that the two dimen- 
sional field <5/*'*(m, v) contains information about the full three 
dimensional structure of the ionosphere if a reference field 

* Only after having introduced this transformation, the author became 
aware of a similar transformation in these publications. Whereas the later 
assume a three-dimensional (i.e. non-planar) array, the focus in the current 
paper is deriving the structure of the three-dimensional (non-planar) iono- 
sphere with a planar array. Combining both is left for a future publication. 



J^q(ii) is available; this is closely related to the technique of 
"holography". We can now go one step further and show after 
a little algebra that 



H(u) = f'^u, v) ® 'K(m, v; w). 



(21) 



being a two-dimensional convolution, with w as control pa- 
rameter. The function 



"Hiu, v; w) = In 



f 

Jo 



e-^'''""^^^Jo(s-U2B)sds, 



(22) 



with M2D = Vm^ + v^ is a Hankel transform of e 2mnVi ^r, s, 
Putting this all together, we find the final result 



m)f = r-' 



6P\u, v) 



/W(m,v) ® 'H{u,v;w) 



(23) 



Hence, the three dimensional power spectrum of the iono- 
sphere or its autocorrelation function can be reconstructed 
from the ratio between the Fourier transform of the rescaled 
scattered intensity field and the Fourier transform of the inci- 
dent radiation field convolved with a Hankel function. 

So what does Eqn.EJ\) mean? Looking carefully at the 
equation, we see that "Hiu, v; w) is the Fourier transform of 
a unit phasor with a phase that is determined by the dis- 
tance (in wavelength) from a half-sphere of radius w to the uv 
plane along a line perpendicular to the latter In other words, 
■?/(«, v; w) act as the Fourier transform of a complex optical 
transfer function in the pupil plane (see e.g. lGoodmanlll98 5). 
which in this case is the full sky and not (as usual) the interfer- 
ometer plane itself. Equivalently, "Hiu, v; w) acts as a complex 
point spread function in uv space (i.e. convolving P'\u,v)). 
Turning this around, a convolution in uv space is a multiplica- 
tion of the sky intensity with the reciprocal of the convolution 
kernel. Hence, 'K(m, v; w) causes a complex beam of unit am- 
plitude on the sky equal to the phasor in Eqn(l22]l. This exactly 
extracts the information from 57*'*(m, v) on a particular w-slice 
cut through the three dimensional autocorrelation function of 
the ionosphere. Hence by simply multiplying the sky model 
I^'\s„, Sv) with the complex phasor in Eqn.(l22li one obtains a 
complex sky-intensity cube. Fourier transforming this back 
slice by slice provides the denominator in Eqn.(l23]l. 

To illustrate this further, imagine that the sky contains only 
one point source of unit flux-density at (so,u, so,v)- In that case. 



nmsr) 



6J^'\u, v) 



-2m(.!o.„u+So.,,i'+ir -0~"^^Io+^' 



(24) 



Substituting Eqn.(fT6]l back into this equation shows that the 
left and right-hand sight of the equation are identical as re- 
quired. By bringing the denominator to the left side of the 
equation, we find for a point source 



|a)(s-S())f =5^-'[w<'>(M,y)], 



(25) 



This is the inversion of Eqn.lfTTTi for a single point source. 
Hence every point source in the sky probes the ionospheric 
power spectrum on an Ewald sphere of reflection. The full 
sky probes the power spectrum of the ionosphere, as encoded 
in Eqn.(l23Tl. on the surface of many offset Ewald spheres of 
reflection, each of unit radius. The function 'H endows each 
point in the sky with a complex phase. In correspondence with 
the Fourier shift theorem (i.e. a phase shift corresponds to a 
spatial shift) this point is then moved in to the third dimension 



Ionospheric Power-Spectrum Tomography in Radio Interferometry 



of power-spectrum which corresponds to the w-direction of 
the ionospher^ 

5.3. The effects of the beam, uv sampling 

Realistic interferometers neither observe the full sky, nor 
sample the electric field in the full plane of the array. If we 
define the complex electric field beam pattern as B(si„ 5,,) and 
the sampling function of the electric field on the plane of the 
interferometer (i.e. usually delta functions at the positions of 
the antennae) as P(u, v) and its Fourier transform as P(su, s,.), 
then Eqn.(l23ll for incoherent emitters can be rewritten as: 

6P^ = [T(|C)|2) X :r(|P|2) X (7« ® -H)] ® ^(^)' (26) 

where we suppressed the explicitly dependence of the func- 
tions on the {su,s^,,Sh) or {u,v,w). The equation holds for 
each given w. The usual two dimensional array point spread 
function (i.e. dirty beam) is denoted by |Pp and its Fourier 
transform is the visibility sampling function in the uv plane. 
Similarly, |B|- is the usual two dimensional antenna beam pat- 
tern. This equation shows clearly that the visibilities of the 
incident intensity field are both multiplied with the Fourier 
transforms of power-spectrum of the ionosphere and that of 
the uv sampling function, respectively. The resulting field 
is then convolved with the Fourier transform of the rescaled 
beam pattern, which to first order (if the beam is much smaller 
than the full sky) is the aperture gain pattern of the antennae. 
To derive |Op) from this equation thus requires deconvolu- 
tion of the observed visibilities of the scattered intensity field. 
This could be difficult in principle were it not for the fact that 
the beam size in general is much larger than the extent of the 
ISP and the dirty beam. In that case the scale over which 
;F(|3)p) X ;F(|Pp) varies is very large compared to the convo- 
lution kernel and the latter can be neglected for modes much 
smaller than the field of view. Hence in practice it is expected 
that the deconvolution is not really required to obtain an ac- 
curate evaluation of the power-spectrum (i.e. |Op). 

6. CONNECTION TO THE PHASE SCREEN 

Whereas in this paper we started our discussion from the 
Born approximation and derived the scattered field, this is 
only a valid approach in the weak scattering scattering or 
weak scintillating regime. Despite this approximation, it is 
remarkably accurate up to intensity fluctuations very close to 
unity. However, t here are other w ays to solve for the scattered 
electric field (e.g. lKravtsovlll992l for a detailed review) in the 
weak scattering regime that are closely connected to the for- 
malism as introduced above and to the phase-screen approach. 
Instead of looking at every point of the medium as a source 
of a single scattering (as is done in the Born approximation), 
one can also assume that the medium does not modify the am- 
phtude of the incident wave to first order and that light rays 

' The phase in the phasor changes by 2n when ws6s « 1 . If the maximum 
of s is set by the field of view and Ss ~ /i/i/max by the resolution of the array, 
the heigh w to which the ionosphere can be probed is Wniax ~ \l(s6s) above 
which phases wrap around multiple times inside a resolution element in the 
image (i.e. information is being lost). Taking a typical field of view of ~5 
degrees and resolution of 3 arcmin (for an array of 2-3 km diameter) leads to 
a maximum high of Wmax'l ~ 25,000 m for /( = 2 m. Structure at ten times 
that height, where the ionosphere typically is densest, is encoded in the image 
on scales ten times below the resolution limit. Hence, baselines well beyond 
20-30 km are required to fully extract the three dimensional structure of the 
ionosphere and its effect on the scattering in the image. One also finds that 
the scale over which the ISP will change is .s ~ Umax/i^'max- This is typically 
tens of degrees for the ionosphere and arc-minutes in the optical. 



travel on a straight line through the medium. In that case, 
since the medium does not absorb or amplify the wave, only 
a phase shift occurs between different points in the medium 
when a plane wave enters the ionosphere. 

We can again describe this in terms of the refractive index or 
electron density as follows. A phase shift of 6tf/ ~ k J 6n(r) ds 
is introduced between a wave traveling in a medium with re- 
fractive index 1 + 6n(r) and unity, respectively, where the in- 
tegral is carried out along a straigh t line through the medium 
(see e.g. 'Luto mirski & Yuralll971h . A plane wave of a source 
with unit amplitude entering medium exits as a wave with a 
"wrinkled" phas e-front due to varying refractive indices along 
different paths (lCronynlll972b . The auto-correlation of this 
"phase-screen" is the function with which the visibilities of 
that source are multiplied in the Mv-plane and its Fourier trans- 
form is the ISP as we discussed it before. This allows us to 
directly connect the previous analysis with that of the phase- 
screen through the so-called Radon transform, which is re- 
lated to the Fourier projection-slice theo rem. 

The Fourier projection-slice theorem (lBracewelllll98 6'), in 
our context, states that the two-dimen sional Fourier transform 
of the Radon projection (Deans 1983), along straight lines, of 
a three-dimensional medium equals a two-dimensional slice 
(perpendicular to the projection direction) through the three- 
dimensional Fourier transform of the object. We now note that 
the phase-screen, up to a constant, is in fact a two-dimensional 
projection of the refractive index 6n{r) in three dimensions. 
Hence the Fourier transform of the phase-screen is simply a 
two-dimensional slice through the three-dimensional Fourier 
transform of the refractive index 6n(r) of the medium. The au- 
tocorrelation of the phase-screen as measured in the Mv-plane 
is then a slice through the three dimensional power spectrum 
of the refractive index, hence that of the electron density dis- 
tribution. We immediately see the connection to the discus- 
sion in Section 5 and how this connects to Eqn. lfTTI) . Since 
point-sources in different directions project the three dimen- 
sional power-spectrum differently on the Mv-plane, they probe 
different slices through the power-spectrum. Disregarding the 
geometric curvature terms, this slice is the tangent plane to the 
Ewald sphere of influence at the point s - So which for small 
vectorial differences are all slices through s - So = (0, 0, 0), 
which is exactly the Fourier projection-slice theorem. 

Observing over a wide field of view allows one to 
build of a three-dimensional electron density power-spectrum 
from these different slices. Obviously the measured auto- 
correlation of the phases are the result of many point sources 
and need to be disentangled. This was discussed in Section 5 
in detail and is identical in the current situation. Thence, the 
phase-screen approach extended to three dimensions is com- 
pletely identical to the approach taken in this paper 

7. RESULTS AND CONCLUSIONS 

A tomographic method has been introduced that allows 
to quantify the three-dimensional power-spectrum of the 
ionospheric electron-density fluctuations based on radio- 
interferometric observations by a two-dimensional planar ar- 
ray. The goal has been to provide a more complete and phys- 
ically intuitive description of the eff'ects of the full three di- 
mensional ionosphere over a wide field of view on radio- 
interferometric images, without any approximations about ei- 
ther a small field of view and/or a thin ionospheric slab (i.e. 
phase screens). Neither is expected to be sufficiently accu- 
rate in upcoming high dynamic range observations with low 
frequency arrays. The description is valid to first-order Born 



8 



Koopmans 



approximation, which holds well for frequencies well above 
the plasma frequency of the ionosphere. Second order cor- 
rections are typically several orders of magnitude below first 
order corrections at frequencies > 100 MHz. However, we 
stress that the method is only valid for weak scattering and 
that higher order corrections are needed if the refractive index 
approaches unity. The main results and conclusions are: 

1. When modeling the ionosphere, it should not be the 
ionospheric electron density distribution that is is pri- 
mary structure to model, but its autocorrelation func- 
tion or equivalent its power-spectrum. Any informa- 
tion about the phase structure of the ionospheric elec- 
tron density distribution is lost in the cross-correlation 
of the electric field when obtaining the visibilities. 

2. A three-dimensional ionosphere causes a spatially vary- 
ing convolution of the sky (at second order level), 
whereas a two-dimensional phase-screen or a thin iono- 
sphere results in a spatially invariant convolution. Iono- 
spheric structure in the w-direction causes, to lowest or- 
der, radial stretching or squeezing of the ionospheric 
scattering point spread function but leaves its tangen- 
tial structure invariant. Correcting for the thickness of 
the ionosphere can thus be reduced (to lowest order) to 
determining a single number, i.e. the level of stretching 
or squeezing. 

3. Residual speckle, which can not be corrected in short 
time integrations, causes a diffuse intensity halo around 
bright sources beyond a certain distance from the 
source. Whereas longer integrations probe more of the 
total scattered flux of the ISP, because of its very steep 
intensity decline away from the source, these longer 
integrations also cause smearing of the instantaneous 
speckle pattern due to variations of the ionosphere. This 
halo ("seeing") and related speckle noise might there- 
fore pose a fundamental limitation on the ability to 
reach the thermal noise level in interferometers at very 
low frequencies after long total integrations. 

4. Long baselines substantially help in correcting for the 
effects of the largest wave modes of the ionosphere, as 
seen inside the field of view of typical interferometers, 
as well as for its three dimensional structure in the in- 
direction. The reason is that these large scale modes 
cause a small-scale diffraction pattern and thus show 
up in the image on small angular scales that can only 
be resolved on long baselines (i.e. of a size equal or 
larger than the projected scale of the wave mode). 

5. An exact mathematical expression is derived that pro- 
vides the power-spectrum of the ionospheric electron- 
density fluctuations from a rescaled scattered intensity 
field and an incident intensity field convolved with a 
complex unit phasor defined on the full sky pupil plane. 
This is related to a holographic principle. In the limit 
of a small field of view, the method reduces to the usual 
thin-phase screen approximation. It is also shown, 
through the application of a Radon projection and the 
Fourier projection-slice theorem that the extension of 
the phase-screen approach to three dimensions is iden- 
tical to the introduced tomographic method. 



Whereas in this paper no direct implementation of an al- 
gorithm is given on how to calibrate the ionosphere, a more 
physically intuitive picture was presented that goes beyond 
the single (or multiple) phase screen models of the ionosphere 
and is valid in the presence of a full three dimensional iono- 
sphere and full-sky field of view. A mathematical expres- 
sion was presented that shows that the ionosphere causes a 
spatially varying point spread function (at the second order 
level) over the field of view, determined by the instantaneous 
three-dimensional power spectrum of the ionospheric elec- 
tron density fluctuations. Subsequently this expression was 
inverted, showing that through a "holographic" principle one 
can extract information about the three-dimensional structure 
of the ionosphere from only electric field measurements on 
a two dimensional plane. In forthcoming work it remains 
to be investigated whether these results can be implemented 
in a practical algorithm in the context of the measurement 
equation (see e.g. Matejek & Morales 2009) in order to cal- 
ibrate visibilities without assumptions about thin ionospheric 
phase screens, and by using instead the full three dimen- 
sional power-spectrum of the ionosphere (i.e. "tomographic 
self-calibration"). 
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